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A spin-charge unifying description for the Hubbard model based on the time dependent local gauge transfor- 
mations is developed. The collective variables for charge and spin are isolated in the form of the space-time 
fluctuating U(l) phase field and rotating spin quantization axis governed by the SU(2) symmetry, respec- 
tively. As a result interacting electrons appear as a composite objects consisting of bare fermions with 
attached U(l) and SU(2) gauge fields. We elaborate on the microscopic origins of the effective action with 
the Coulomb interaction that contain topological theta terms. Furthermore, we unravel the link between non- 
trivial multiply-connected topological structure of the U(2)=U(1)®SU(2) configurational space for gauge 
fields and the instanton contribution to the statistical sum. 
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1 Introduction 

Among the electronic Hamiltonians relevant for interacting systems the Hubbard model (U El is con- 
sidered as one that contains the essential ingredients for understanding the physics of correlated electrons. 
In a many-body system the relevant physics is encoded in the symmetries and in the Hubbard model they 
are represented by the charge U(l) gauge and spin rotational SU(2) groups relevant for the occurrence 
of the superconducting and magnetic orderings. Interestingly, the resulting U(2)=U(1)®SU(2) symme- 
try group has a nontrivial multiply connected topological structure. For multiply connected configuration 
spaces novel features can arise as for example in the Aharonov-Bohm effect [4] governed by the multiply 
connected U(l) group manifold. Since the homotopy class 7Ti[U(2)]= Z also forms a set of integer wind- 
ing numbers the topological structure of the configuration space is nontrivial, ambiguities may arise when 
attempts are made to specify a value for the phase of a wave function for the whole configuration space 
0. Thus the problem we are facing is that of many -body quantum mechanics on a multiply connected 
configuration space. In the present work we develop a spin-charge unifying description for interacting 
electrons given by the Hubbard model. It is based on the time dependent local U(2) gauge transforma- 
tions to disentangle the Coulomb interaction. The collective variables for charge and spin are isolated in a 
form of the space-time fluctuating U(l) phase field and the rotating spin quantization axis governed by the 
SU(2) symmetry, respectively. As a result interacting electrons appear as composite objects consisting of 
bare fermions with attached U(l) and SU(2) gauge fields and due nontrivial topology of the U(2) group an 
instanton contribution to the statistical sum emerges. 
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2 Hubbard Hamiltonian 

Our starting point is the purely fermionic Hubbard Hamiltonian H = Ht + Ti-u'- 

U = - t E [4(r)c a (r')+b-.c.] + C/^n T (r)n ; (r). (1) 

(rr'),a r 

Here, (r, r') runs over the nearest-neighbor (n.n.) sites, t is the hopping amplitude, U stands for the 
Coulomb repulsion, while the operator c' a (r) creates an electron with spin a =f, J, at the lattice site r, 
where n a (r) = cJ l (r)c Q ,(r). Usually, working in the grand canonical ensemble a term is added to 7i 
in Eq.(fl]ito control the average number of electrons, TL — > H — nJ2 r n ( r ) w ' trl A* being the chemical 
potential and n(r) = u-j-(r) + ^j.(r) the number operator. It is customary to introduce Grassmann fields, 
c Q (rr) depending on the "imaginary time" < r < (3 = 1/kgT, (with T being the temperature) that 
satisfy the anti-periodic condition c a (rr) = — c a (rr + f3), to write the path integral for the statistical sum 
Z = J [VcDc] e^ s ^ c ^ with the fermionic action 

rP 

S[c, c] = Sb[c, c] + / (1tH[c,c\, (2) 
Jo 

that contains the fermionic Berry term (6): Sb [c, c] = J2ra fa dTc a {YT)d T c a {vT) that will play an im- 
portant role in our considerations. 



3 Spin-charge U(2) reference frames 

For strongly correlated system it is crucial to construct a covariant formulation of the theory which naturally 
preserves the spin-rotational symmetry present in the Hubbard Hamiltonian. In order to maintain spin 
rotational invariance, one should consider the spin-quantization axis to be a priori arbitrary and integrate 
over all possible directions in the partition function. For this purpose the density-density product in Eq.(Q~|i 
we write, following Ref.flT), in a spin-rotational invariant way: 



where S' a (rr) = \ ^2 aa i c l( rr ) | 5'aa' c «'( rr ) denotes the vector spin operator (a = x, y, z) with a a being 
the Pauli matrices. The unit vector r2(rr) = [sini^rr) cos <p(rr), sin i?(rr) sin f{rr), cos i9(rr)] written 
in terms of polar angles labels varying in space-time spin quantization axis. The spin-rotation invariance is 
made explicit by performing the angular integration over J7(rr) at each site and time. By decoupling spin 
and charge density terms in Eq.© using auxiliary fields q{vt) and iV(rr) respectively, we write down the 
partition function in the form [9 | 

Z = J[Dfi\ J[VVVg] J [VcDc]e- s[n < v ' e > 5 < c] . (4) 
where [Dfl] = Y[ rT S '" f>tyTTk ^ d ^ rTk ^ d[p ^ rTk ^ is the spin-angular integration measure. The effective action 



lrT fc Ait 

reads: 



S[n,V,g,c,c] = / dT 

r- JO 



2 {yt) V 2 (yt) 

^j- 1 + — fj- 1 + iV{Yr)n(YT) + 2q(yt)VL{yt) ■ S(rr) 



r P 

S B [c lC }+ I dTU t {c,c\. (5) 



Simple Hartree-Fock theory won't work for a Hubbard model in which U is the largest energy in the prob- 
lem. One has to isolate strongly fluctuating modes generated by the Hubbard term according to the charge 
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U(l) and spin SU(2) symmetries. To this end we write the fluctuating "imaginary chemical potential" 
iV(rr) as a sum of a static Vb(r) and periodic function V(tt) = Vb(r) + V(rr) using Fourier series 
V{yt) = I J2n=iiV( ruJ n)e iuJnT + c.c] with w B = 27m//3 (n = 0, ±1, ±2) being the (Bose) Matsubara 



frequencies. Now, we introduce the U(l) phase field 
_ d<j>(rr) 1 d 



(rr) 



r) -f e 



!>(rr) via the Faraday-type relation 1 8 ] 
<« rr ) = V{vt). 



Furthermore, by performing the local gauge transformation to the new fermionic variables / Q (rr): 



c Q (rr) 




_ c Q (rr) 





z(rr) 
z(rr) 



/a(rr) 
/a(rr) 



(6) 



(7) 



where the unimodular parameter z(rr)| = 1 satisfies z(rr) = e 



we remove the imaginary term 



i J drV (rr)n(rr) for all the Fourier modes of the V^(rr) field, except for the zero frequency. Subsequent 
SU(2) transformation from f a (rr) to h a (rr) operators, 



" /t(rr) " 




' C T (rr) 


-Cl(rr) " 




' fci(rr) " 






. Ci(rr) 


Ct(rr) 







R(rr) = 



C T (rr) 
Ci(rr) 



-Ci(rr) 
Ct(rr) 



(8) 



with the constraint |(f(rr)| 2 + |£|(rr)| 2 = 1 takes away the rotational dependence on l~i(rr) in the 
spin sector. This is done by means of the Hopf map R(rr)tj z Rt(rr) = a ■ fi(rr) that is based on the 
enlargement from two-sphere Si to the three-sphere S3 ~ SU(2). The unimodular constraint can be 
resolved by using the parametrization 



Cl(rr) 



= g-§[v( rr )+x(r 



")], 



e i[ V (rr)- X (rr)] sin 



i?(rr) 
2 

i?(rr) 



(9) 



with the Euler angular variables v?( rT )> $(rr) and x( rT X respectively. Here, the extra variable x( rr ) 
represents the U(l) gauge freedom of the theory as a consequence of S% — ► S3 mapping. One can sum- 
marize Eqs (Q and (O by the single joint gauge transformation exhibiting electron operator factorization 
c a (rr) = ^2 a i U aa i (YT)h a i (rr), where U{yt) = z(rr)R(rr) is a U(2) matrix which rotates the spin- 
quantization axis at site r and time r. This reflects the composite nature of the interacting electron formed 
frombosonic spinorial and charge degrees of freedom given by R aa ' (rr) and z(rr), respectively as well as 
remaining fermionic part h a (rr). Accordingly, the integration measure over the group manifold becomes 

2t 0(r/3)=0 o (r)+27rm(r) f2(r/3)=f2 

J [V^VSl] = J2 II / d ^ j dfl °^ J V4>{vt) J PO(rr) (10) 

where J dfl ■ ■ ■ = j- f~ sin Odd dip . . . and [2?S~2(rr)] = ]~| fe dfl(rrk). Here, m E Z labels equiva- 
lence classes of homotopically connected paths |[T0l . 



4 Effective phase-angular action 

Anticipating that spatial and temporal fluctuations of the fields Vb(r) and g(rr) will be energetically pe- 
nalised, since they are gaped and decouple from the angular and phase variables. Therefore, in order to 
make further progress towards we next subject the functional to a saddle point analysis. The expectation 
value of the static (zero frequency) part of the fluctuating electrochemical potential Vb (r) we calculate by 
the saddle point method to give Voir) — i (/i — |fi),) = ip where ft is the chemical potential with a 
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Hartree shift originating from the saddle-point value of the static variable Vb(r) with nh = nh\ + nui and 
nha = (/i a (rT)/i a (rr)). Similarly in the magnetic sector 

{ - l lk 

where A c = U(S z (rT) sets the magnitude for the Mott-charge gap. The two choices delineated in Eq.dTTl 
correspond to the saddle point of the "antifferomagnetic" (with staggering A c ) or "ferromagnetic type". 
Note that the notion ferromagnetic (antifferomagnetic) here does not mean an actual long-range ordering - 
for this the angular spin-quantization variables have to be ordered as well. In the new variables the action in 
Eq.© assumes the form S [fi, g, h, h] = Sb [h, h] + J Q drrin,^ [p, h, h] + So [<p] + 2 ^ r J* dr A(rr) ■ 
S h (rr), where S h (rr) = ± ^a 7 ^(r'O'^M 1 '''")- Furthermore, 



. Jo 



(f 2 rr 12//;. , 
— - — - H -6(tt) 



(12) 



stands for the kinetic and Berry term of the U(l) phase field in the charge sector. The SU(2) gauge trans- 
formation in Eq.© and the fermionic Berry term generate SU(2) potentials given by R^(rr)9 r R(rT) = 

Rt + #m + *&) R = • A ( r ^)' where 

A x (vt) = -i?(rr) sin x(rr) — -(p(rr) sin 0(rr) cos x(rr) 

i . i 
A y (rr) — -i?(rr) cos x( rr ) + ~^ l f( rr ) s i n 0(rr) sin x(vr) 

i i 

A z (rr) = -<p(rr) cos^(rr) + -x(rr). (13) 



are the SU(2) gauge potentials. The fermionic sector, in turn, is governed by the effective Hamiltonian 



rin,4> = Y] e( r T)[M rr )M rr ) - ^(rr)^(rr)] 



-'EE [W'frrMr'rJl^^H^tr'rj^E^t^^H. ( 14 ) 

(r.r') ck7 ro 

The result of the gauge transformations is that we have managed to cast the strongly correlated problem 
into a system of mutually non-interacting fermions, submerged in the bath of strongly fluctuating U(l) and 
SU(2) fields whose dynamics is governed by the energy scale set by the Coulomb interaction U coupled to 
fermions via hopping term and with the Zeeman-type contribution with the massive field g(rr). 

In analogy to the charge U(l) field the SU(2) spin system exhibit emergent dynamics. Integration of 
fermions will generate the kinetic term for the SU(2) rotors 



1 ^ 

S [Q] = x 4 yW drdr' V A a (rr)A fc (r'r') 

2 Jo ab 



£<S£(rr)S£(rV)) - 2]T f drA{vr) ■ <S fc (rV)> (15) 
/ i Jo 



with 



(S a h (rr)S b h (r'r')) = -\x2S ab j-, (S* h (rr)) = ^ (16) 

In the AF phase, at the half-filling, it assumes the staggered form g(rr) = A c (— l) r with A c being the 
charge gap A c ~ U/2 for U/t ;» 1, with — \J e\ + A;?, Q = (n, it) and e k = e k — p,. At zero 
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1 A 

temperature liniT^o ■g- — ^gS" so l ^ at ' m me limit U/t ^> 1 one obtains £ s ~ 2A C = {/. However, 
a nonzero value of A c does not imply the existence of anfiferromagnetic long-range order. For this the 
angular degrees of freedom fi(rr) have also to be ordered, whose low-lying excitations are in the form of 
spin waves. Therefore the kinetic term in the spin sector becomes <Sq [SI] = -j-Er/o <iTA(rr) • A(rr) 
so that the total kinetic energy S[Sl, <j>] — Sq[4>] + Sq[SI] for U(l) and SU(2) rotors becomes 




12/i- 1 



+ 2<y9(rr)x(rr) cos$(rr)] + ^ [0(rr) cosz9(rr) + x(rr)] | . (17) 

The distinctive feature of Eq.(TT7b is th presence of the geometric Berry contributions that signify topolog- 
ical features of the underlying field theory ifTTI . 



5 Summary 

In this work we have observed that the important symmetries of the Hubbard model given by by the charge 
U(l) gauge and spin rotational SU(2) groups imply that the quantum dynamics is governed non-trivially 
by the multiply connected U(2)=U(l)<g>SU(2) group manifold. Following this route we have developed 
a spin-charge unifying description for interacting electrons given by the Hubbard model. It is based on 
the time dependent local gauge transformations to disentangle the Coulomb interaction. The collective 
variables for charge and spin are isolated in a form of the space-time fluctuating phase field and the ro- 
tating spin quantization axis. As a result interacting electrons appear as a composite objects consisting of 
bare fermions with attached gauge fields. Our results form a general framework tailored to describe eg. 
both magnetic and superconducting order which is rooted in microscopic of strongly correlated electrons 
and their basic symmetries and interactions lfT2l . In particular, it will enable to identify the degrees of 
freedom of the hierarchy of low-energy effective theories, and the way they are associated with the global 
symmetries. 
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